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ABSTRACT 

We investigate the physical properties of dense cores formed in turbulent, 
magnetized, parsec-scale clumps of molecular clouds, using three-dimensional nu- 
merical simulations that include protostellar outflow feedback. The dense cores 
are identified in the simulated density data cube through a clumpfind algorithm. 
We find that the core velocity dispersion does not show any clear dependence 
on the core size, in contrast to Larson's linewidth-size relation, but consistent 
with recent observations. In the absence of a magnetic field, the majority of the 
cores have supersonic velocity dispersions. A moderately-strong magnetic field 
reduces the dispersion to a subsonic or at most transonic value typically. Most of 
the cores are out of virial equilibrium, with the external pressure dominating the 
self-gravity. The implication is that the core evolution is largely controlled by 
the outflow-driven turbulence. Even an initially-weak magnetic field can retard 
star formation significantly, because the field is amplified by the outflow-driven 
turbulence to an equipartition strength, with the distorted field component dom- 
inating the uniform one. In contrast, for a moderately-strong field, the uniform 
component remains dominant. Such a difference in the magnetic structure is 
evident in our simulated polarization maps of dust thermal emission; it provides 
a handle on the field strength. Recent polarization measurements show that the 
field lines in cluster-forming clumps are spatially well-ordered. It is indicative 
of a moderately-strong, dynamically important, field which, in combination with 
outflow feedback, can keep the rate of star formation in embedded clusters at the 
observationally-inferred, relatively- slow rate of several percent per free-fall time. 

Subject headings: ISM: clouds — ISM: magnetic fields — MHD — polarization — 
stars: formation — turbulence 
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Introduction 



Millimeter and submillimeter observations of dense cores in nearby parsec-scale 
cluster-forming clumps have shown that the core mass functi o n (CM F) resembles the 



MotteetaL 



(Il998[ ) performed 1.2 mm 



stellar initial mass function (IMF). For example, 
continuum observations toward L1688 and identified 57 starless dense cores whose mass 
function {dN/dM ~ M~^-^) is consistent with the Salpeter power-law IMF at the high 
mass end {dN/dM ~ M~^-^^). The CMF also has a break at about 0.3 Mq, below which it 



flattens to dN/dM ~ M . Other authors have found simi 



Johnstone et al 



2000 : 



Stanke et al 



2006 



Maruta et al 



ar CMFs in the region (e.g. 



20101). The sha pe of the CMF i s 



broad ly consistent with the IMF of Class II YSOs in the same region (iBontenips et al. 



Testi fc SergentI ( 



19981 ) identified, 



20011 ). Another example is the Serpens Cloud Core, where 
through 3 mm continuum interferometric observations, 32 dense cores whose mass function 
can be fitted by a power-law of dN/dM ~ M~^-^, again consistent with the Salpeter IMF. 
These and other observations suggest that the identified dense cores may be the direct 
progenitors of individual stars and the bulk of the stellar IMF may be at least partly 
determined by cloud fragmentation in the parsec-scale dense clumps. Thus, understanding 
the formation process of dense cores is a key step towards a full understanding of how stars 
form. 

Both supersonic turbulence and magnetic fields are expected to play a role in clustered 
star formation. The relative importance of the two is still under debated, however. One 
school of thought is that the core (and thus star) formation in a cluster-forming clump is 
mainly regulated by supersonic turbulence. In this scenario, star formation is completed 



rapid 



2000 



, before the initial supersonic turbu 



Hartmann et al. 



picture, although 



2001 



Klessen et al. 



Padoan fc Nordlund 



ence has decayed significantly (e.g., lElmegreen 



20001 ). Magnetic fields are less important in this 



( I2OO2I ) argued that a weak magnetic field is still 
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needed in order to produce a core mass spectrum that resembles the stellar IMF. The weak 
magnetic field is expected to be strongly tangled by the supersonic turbulence. Ordered 
field structures, if present at all, are typically attributed to large-scale compression and 
are expected to b e parallel to the dens e filaments that are also produced by the same 



compression (e.g.. 



Pelkonen et al. 



20071 ) 



The second school of thought envisions a magnetic field that is more dynamically 
important; it prevents stars from forming too rapidly. Because supersonic turbulence 
dissipates rapidly, it must be somehow replenished. In a cluster-forming clum p, protostellar 



outfiow feedback can p 



2009 



Arce et al. 



2010 



ay an important role in tu rbulence regeneration ( e.g. 



Nakamura et al. 



2011aU bl). iLi fc Nakamural (120061 ) demonstrated 



Maury et al. 



that protostellar outfiows can resupply the supersonic turbulence, keep ing the clumps 



near a quasi-virial equilibrium state for a relatively long time (see also 



Nakamura &: Li 



2007 



Matzner 



2007 



Carroll et al. 



20091 ). Because the matter moves preferentially 



along the dynamically important magnetic field, the field is expected to be more or less 
perpendicular to the dense filaments that are created by either turbulence compression or 
self-gravity. A goal of this paper is to quantify the effects of the magnetic field on the 
clustered star formation in general, and the properties of dense cores in particular, through 
3D MHD simulations that include outfiow feedback for turbulence replenishment. 

The rest of the paper is organized as follows. In Section [2], we describe the numerical 
method and simulation setup. In Section [3l we apply a version of the so-called "clumpfind" 
algorithm to identify dense cores in a parent dense clump and derive their properties. We 
find that the internal motions of the cores are sensitive to the magnetic field strength of the 
parent clump and that the external surface pressures play a dominant role in core formation 
in outfiow-driven turbulence. In Section HI we compare the properties of the simulated cores 
with those observed in a couple of nearby cluster forming regions. Our main conclusions 
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are summarized in Section |5l 



Model Formulation 



The initial and bou n dary conditions of the simulations are essentially the same as 



those of iNakamura fc Li 



(120071 ). We consider a centrally-condensed molecular gas clump 
inside a cubic simulation box of L = 1.5 pc on each side, with an initial density profile 
of p{r) = Po[l + [r/rcY]^^ for r < L/2 (where r^ = L/6 is the radius of the central 
plateau region) and p{r) = O.lpo for r > L/2. Here, the central density po is given by 
Po = 4.68 X lO^^^nHa.o g cm~^, with nH2,o being the central number density of molecular 
hydrogen, assuming 1 He for every 10 H atoms. Periodic boundary condition is applied to 
each side of the box. We adopt a central H2 density of n^/j.o = 2.69 x 10'^(T/20 K)(1.5 pc/L)^ 
cm~^, corresponding to a central free-fall time tg^ = 0.19 Myr and a central Jeans length of 
Lj = {Ticl/GpoY/^ - 0.17(T/20K)V2(riH,,o/2.69 x 10^cm-3)-i/2 pc. it yields a total clump 
mass of Mtot = 939M0. The average clump density is nnj = 4.0 x 10'^ cm~^, corresponding 
to a global free-fall time tg^ci = 0.49 Myr. An isothermal equation of state is assumed, with 
a sound speed of Cg = 0.266 km s^^ for gas temperature T = 20 K. 

Our choice of the clump mass, density and gas temperature is motivated by observations 

orming clump, the p Ophiuchi main cloud. From ^'^CO 



of the nearest pc-scale cl uster 

( J = 1 — 0) observations, 

L1688 region, adius ted for the parallax distance of 125 pc (e.g.. 



LorenI (119891 ) estimated a total gas mas s of 865 Mf^ for the w hole 



Lombardi et al. 



2008 



Loinard et al. 



20081 ). The gas temperature in the region has some spatial variation, but 



the average appears close to 20 K. O ur clump paramet ers are also consistent with those of 



other nearby cluster-forming clumps (JRidge et al. 



2OO3I ) and the i nfrared dark c 



are thought to be future sites of cluster formation in GMCs (e.g. 



Butler fc Tan 



ouds which 



20091). 
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At the beginning of the simulation, we impose on the clump a uniform magnetic field 
along the x-axis. The field strength is specified by the plasma /3, the ratio of thermal to 
magnetic pressures at the clump center, through Bq = 4 7/3~^/^(T/20K)^/^(nH„ n/2.69 x 
10^cm-3)i/2 ^G. In units of the critical value l/{27iG^/'^) jNakano fc NakamuralllQTsJ ). the 
mass-to-fiux ratio in the central flux tube is given by Aq — 8.3/3^'^. The mass-to-flux ratio 
for the clump as a whole is A ~ 3.0/3^/^. Although systematic measurements of magnetic 
field strength in cluster-forming clumps are not available, relatively strong magnetic fields 



are inferred in some cases: e.g., 850 /xG for OMC 1 at ra 



Houde et al 



2004' 



( ISugitani et al. 



Nakamura Sz Li 



I 



10^ ^ cm ^ (see 



Crutcher 



1999 



2007f ) and 160 /iG for Serpens at riHa.o ~ 10^~^ cm~^ 
2010l ). In the present paper, we adopt as a representative value /3 = 0.2 (or 
100 fiG at nHa ~ 3 x 10^ cm~^), corresponding to A ^ 1.4 or a dimensionless fiux-to-mass 
ratio r ~ 0.74, so that the clump as a whole, as well as the denser central region, is 
magnetically supercritical. 



Following the standard practice (e.g. 



Mac Low et al 



1998; 



Ostriker et al 



200l|), we 



stir the initial clump at the beginning of the simulation with a turbulent velocity field 
of power spectrum v^ oc k^"^ and rms Mach number A4 = 5. Our initial clump has a 
virial parameter (the ratio of the kinetic energy to the gravitational energy) of Ovir — 0.5. 
The effective Ovir should be close to unity for T ^ 0.74 since the gravitational energy is 
effectively reduced by the magnetic field by a factor of 1 — F^. The turbulence is allowed to 
decay freely, except for feedback from protostellar outflows. 

The evolution of the turbulent, magnetized molecular clump is followed using a three- 
dimensional MHD cod e based on an upwind TVD scheme. The MHD cod e is essentially the 



same as those used in 



Li fc Nakamural (120061 ) and 



Nakamura fc Li 



(120071 1. The ideal MHD 



equations are solved using the code having second-order accuracy in both space and time. 
The equation of state is assumed to be isothermal. To ensure the divergence-free magnetic 



-7- 

field, the divergence cleaning method is adopted. The Poisson equation for gravitational 
potential is solved using the fast Fourier transform. Our simulation has a resolution of 
256^. Although our focus is on core properties, we do inclu de a crude treatment of the 



formed stars and outflows, following 



Nakamura fc Li 



(120071 ): we refer the reader to that 



paper for details. Briefly, when the density in a cell crosses the threshold pth = 400^0, 
corresponding to 10^ cm^'^, we create a Lagrangian particle at the cell center. We extract 
mass from a small region surrounding the cell, and put it on the particle, which moves 
with mass-weighted mean velocity of the extracted mass (see Nakamura & Li 2007 for more 



detail). After creation, the particle is a. 



to the same prescription as 



Wang et al. 



lowed to accrete from the surrounding gas according 



(J2010l ). To mimic the effect of protostellar outflows. 



the particle injects into the ambient gas a momentum that is proportional to the particle 
mass increment AM*. The outflow momentum is scaled with the dimensionless outflow 
parameter / as P = f(AM^/Mf7^)(V ,„/100kin s^^), where V^ is the wind velocity (see the 



Appendix of 



Nakamura &: Li 



20071 ). The fiducial value of / is set to 0.4, consistent with 



Matzner fc McKed (J2000h . Each outflow has a bipolar and spherical component, with a 



momentum ratio of 3:1. The model parameters are summarized in Tabled! 



Numerical Results 



We first concentrate on the dense cores identified in three representative models: model 
Nl (no magnetic field, /3 = oo), Wl (weak magnetic field, (3 = 2), and SI (moderately-strong 
magnetic field, (3 = 0.2). Their dimensionless outflow parameter is set to the fiducial value 
/ = 0.4. 



3.1. Core Identification 



Following 



Nakamura fc Lil (J2008[ ). we identify dense cores from the d ensity data cube 



Williams et al 



fll994l ). We adopt a 



using a variant of the CLUMPFIND algorithm of 
threshold density pth,min = 4po (~ 10^ cm"'^), and a maximum density pth.max = 400 po 
(~ 10^ cm~^, above which a Lagrangian particle is created). Thus, the cores identified 
from our simulation data correspond roughly to the dense cores observed in the dust 
continuum emission and high density molecular tracers such as N2H+ (J = 1 — 0) and 
H^'^CO^ {J = 1 — 0). The density distribution between pth,min and pth,max is divided into 
10 bins equally spaced logarithmically. We include only those spatially resolved cores 
containing more than 50 cells in the analysis, to ensure that their properties are determined 
with reasonable accuracy. We also choose the cores that do not include particles. In 
this sense, our identified cores correspond to starless cores. The minimum core mass and 
radius identified through this procedure are Mc^min = Pth.min x SOAxAyA^; ~ O.OSM© 
and ~ 0.01 pc, respectively. We note that our core identification procedure is not 
exactly the same as those used in observational studies, which are typically based on 
position-position-velocity data cubes of molecular line emission or column density maps 
from dust continuum emission. 

3.2. Spatial Distribution of Dense Gas 

Before discussing the physical properties of the identified cores, we show how the 
magnetic field affects the global density distribution. Figure [U compares the column density 
distributions along the y-axis for the three models with different initial magnetic field 
strengths, at a stage when the star formation efficiency has reached 16%. Since the star 
formation is retarded by the magnetic field, the evolution time tends to be longer at a 
given SFE for a stronger magnetic field. For comparison, the positions of the formed stars 
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and the identified dense cores projected on the x — z plane are indicated by the dots and 
triangles, respectively, in Fig. [TJ The initial magnetic field direction is parallel to the x-axis 
(the abscissa). 

Figure [1] indicates that the global density distribution depends on the initial magnetic 
field strength. For the non-magnetic model (model Nl), the global density distribution 
shows a large-scale filamentary structure that contains many small fragments (Figs. [1^ and 
[T]i). The formation of the large-scale filament is probably because the turbulence energy 
is initially highest on the largest scale and the large-scale compression may be further 
amplified by self-gravity to create the filamentary structure. 

For the weakly-magnetized case (model Wl), the global density distribution does 
not clearly show a large-scale filamentary structure by the stage shown in Figs. [TJj and 
[1^. The dense part appears to contain many smaller filaments or elongated fragments 
that are distributed almost independently of the initial magnetic field direction. This is 
different from the non-magnetic case. The difference may be because the local turbulent 
motions amplify random magnetic field components (see Section H^ . which slow down the 
formation of the large-scale filamentary structure. 

In contrast, in the presence of a moderately-strong magnetic field (model SI), the 
cloud material condenses preferentially along the magnetic field lines into a large-scale 
filamentary structure that is nearly perpendicular to the initial magnetic field direction. 
Several less-dense filaments are also associated with the large-scale filament (hereafter 
the main filament). The less-dense filaments appear to be more or less along the initial 
magnetic field direction and to merge toward the dense parts of the main filament. 
The dense cores are distributed primarily along the main filament. The column density 
dis tribution is re miniscent of the hub-filament structures of star forming regions pointed out 



by 



Myerd ( 120091 ). The observed hub-filament structures may indicate the presence of the 
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moderately-strong, spatially-ordered magnetic field, instead of the weak magnetic field that 
are tangled by supersonic turbulence. We will revisit this issue in Section 14.3^ where we 
compare dust polarization maps derived from the simulation data for models Wl and SI. 



3.3. Core Property 



As shown in 



Li fc Nakamural (120061 ) and iNakamura fc Lil (120071 ) . our model clump has 



reached a quasi-equilibrium state within one global free-fall time. After that, the statistical 
properties of the cores remain largely-unaltered. Therefore, to compare the statistical 
properties of the identified cores, we choose the same stage as that of Fig. [H i.e., the 
stage when 16 % of the total mass has been converted into stars. In this subsection, we 
compare the core properties for the three representative models, models Nl, Wl, and SI. 
The total numbers of identified cores are 160, 105, and 105 cores for models Nl, Wl, and 
SI, respectively. 

3.3.1. Radius, Mass, and Velocity Dispersion 

We now focus on the physical properties of the cores identified in the three representative 
models. To facilitate comparison with observations, we present the core properties in 
dimensional units. Figures Et through |2l show the histograms of the core radius, mass, and 
velocity dispersion of the nonthermal component. The minimum, maximum, mean, and 
median values are summarized in Table [21 The core radius, -Re, is defined as the radius of 
a sphere having the same volume as the core. The distributions of core radius are broadly 
similar in all cases. They range from ~0.01 pc to ~0.05 pc, all peaking around 0.02 pc. The 
maximum and mean core radii are somewhat larger in model SI (the moderately-strong 
field case) than in the other two. The core mass distribution shows a similar trend. In all 
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models, the core mass ranges from ~ 0.07 Mq to ~ AMq. The mean value is somewhat 
larger for model SI than for the other two. The mean core mass and radius are larger in 
model SI probably because the stronger magnetic field can support more mass against 
gravitational collapse. 

The magnetic effect is even more prominent in the distribution of core velocity 
dispersion. For the non- magnetic model (model Nl), the velocity dispersion ranges from 
~ 0.1 km s~^ to 4 km s^^, peaking around 1 km s~^ (corresponding to a Mach number 
of ~ 4). Most of the cores have supersonic internal motions in this case, which appears 
inconsistent with the predominantly subsonic or at most transonic internal motions inferred 
observationally (e.g., Andre et al. 2007; Walsh et al. 2007; Maruta et al. 2010, see also 
Section |4?T|) . Even a relatively weak magnetic field appears to affect the internal motions 
inside the cores significantly. For model Wl, the median velocity dispersion is estimated at 
~0.5 km s~^, corresponding to a Mach number of about 2, although the mean is almost 
the same as that of model Nl. Most of the cores in this case still have supersonic velocity 
dispersions. 

In contrast, for model SI, the majority of the cores have a velocity dispersion close 
to, or less than, ~ 0.3 km s^^, which is comparable to the sound speed (Cg = 0.266 km 
s~^). The velocity dispersion ranges from ~ 0.05 km s~^ to 2 km s~^, peaking around 
the median of 0.3 km s~^. In other words, the strong magnetic field tends to reduce the 
core internal motions significantly. The reduction is probably because the magnetic field 
cushions the converging flows in the supersonic turbulence that are largely responsible for 
the core formation. The magnetic effect is also evident in the virial state of the cores, as we 
show below. 
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3.3.2. Velocity Dispersion-Radius Relation 

The internal motions of the observed dense cores are often measured using the velocity 
dispersions derived from molecular line data. Here, we compare the velocity dispersions 
for the three representative models with different field strengths to study the effects of the 
magnetic field. 

Figure [3] plots the non-thermal components of 3D velocity dispersions of the identified 
cores against the core radii for (a) model Nl, (b) model Wl, and (c) model SI. The 
best-fit power-laws are given by (A^Nx/km s-^) = (2.2 ± 1.9)(i?c/pc)°-^^^°-22 (7^ = 0.023), 
(A\/NT/km s-i) = (3.8 ± 9.0)(i?c/pc)°-^^^°-^^ (7^ = 0.070), and (A^NT/km s'^) = 
(1.1 ± 0.48)(i?c/pc)-°-3^^°-^2 (7^ = 0.29) for models Nl, Wl, and SI, respectively, where 7^ 
denotes the correlation coefficient. It shows that there is no clear correlation between the 
velocity dispersion and core radi us. Our cores do not follow the power-law linewidth-radius 



relation derived by 



Larson 



(jl98ll ) . The same conclusion holds for dense core s formed in even 



more strongly magnetized, magnetically subcritical, clouds (see Figure 18 of 



Nakamura &: Li 



20081 ) . The reason is probably that the dense cores created out of turbulent molecular gas 
are generally not in virial equilibrium (the equilibrium is implied by Larson's laws) and 
that their formation is strongly influenced by external turbulent flows (see below). 

Quantitatively, the velocity dispersions of the cores depend on the initial magnetic field 
strength, as shown already in Figs. Wg through [2}. In the absence of magnetic field (model 
Nl), about 40% of the cores have velocity dispersions larger than 1 km s~^; the majority 
of the cores have supersonic internal motions. As mentioned above, even a relatively weak 
magnetic field can influence the internal motions significantly. In model Wl, although many 
cores still have large velocity dispersions, the median value is less than that in model Nl by 
a factor of ~ 2. In the stronger field case (model SI), the median and mean of the velocity 
dispersions are estimated at 0.31 and 0.46 km s^^. About 40 % of the cores have velocity 
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dispersions smaller than the sound speed of 0.266 km s~^; the majority of the cores have 
subsonic or at most transonic internal motions, which is more consistent with observations. 
We conclude that the magnetic field plays a significant role in reducing the internal motions 
of the cores. 

3.3.3. Virial Parameter 

The virial parameter, the ratio of the virial mass to core mass, is often used as a 
measure of the gravitational boundedness of dense molecular cloud cores, particularly in 
observational studies. In Fig. IU we plot the virial parameters against the core mass for the 
three models. Here, we calculate the virial parameter as 

a . = 5^^iB^- (1) 

where AV^id is the ID FWHM velocity width including thermal contribution and the virial 
mass of each core is obtained by assuming that it is a uniform sphere. The virial mass 
depends on the density distribution. For a centrally condense sphere with p oc r~^, the 
virial parameter is reduced by a factor of 5/3. Note that the effects of the nonspherical 
mass distri bution appear to be smal l as long as the aspect ratios of the cores are not far 



from unity (JBertoldi fc McKee 



1992). 



Figure H] indicates that the relationship between the virial parameter and core mass 
depends on the initial magnetic field strength. In the absence of a magnetic field (model 
Nl), the virial parameters are typically very large, and have a large scatter. The scatter 
in the virial parameter-mass plot comes mainly from the large scatter in the velocity 
dispersion-radius relation (see Fig. [3]). The virial parameters often exceed 10 — 10^ even 
for relatively massive cores, indicating that, for most of the cores, the gravitational energy 
is much smaller than the internal kinetic energy. The virial parameters tend to increase 
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with decreasing core mass, and the lower bound of the virial parameters follows roughly a 



power-law relation of ayir oc Mc 



-2/3 



For the model with a weak magnetic field (model Wl), the virial parameter- mass 
relation is qualitatively similar to that of the non-magnetized model (model Nl). Although 
the fraction of cores with large virial parameters is somewhat small for the weak magnetic 
field case, the scatter in the viral parameter-mass relation remains large. About 30 % of the 
cores still have virial parameters larger than 10. 

On the other hand, for the model with a moderately-strong magnetic field (model 
SI), the identified cores tend to follow a single power-law relation of a^a oc Mc , with 
a rather small scatter. Furthermore, most of the cores (about 9 % of the cores) have 
virial parameters smaller than 10 . This trend is consis tent with the initially-magnetically 



subcritical case (see Figure 17 of 



th e scaling found by 



to 



Nakamura fc Li 



20081 ). It is also in good agreement with 



Bertoldi fc McKed (Il992[ ) for nearby giant molecular clouds. According 



Bertoldi fc McKed ( 1l992h . an object confined by the ambient pressure has a virial ratio 



that is p roportional to {MjMj) ^ /^, where Mj is the Jeans mass defined by Equation 



(2.13) of 



Bertoldi &: McKed ( 1l992l ). The fact that our cores follow the same relation implies 



that the ambient pressure plays an important role in confining our cores as well. In the 
next subsection, we perform a detailed virial analysis of the cores to quantitatively clarify 
the role of the surface pressure term in core dynamics. 

3.3.4- Vinal Analysis 



The virial theorem is useful for analyzing the core dynamical properties. Here, we 
perform a detailed virial analysis using the virial equation in Eulerian coordinates (e.g.. 
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McKee fc Zweibe]|ll992l : iTillev fc Pudritzl 120071 : bib et all 120071 ) 

S = -! + -— / pr^v ■dS = U + K + M + W + S + F 
2 2 dt J (j 



(2) 



where / = jypr'^dV, U = SjyPdV, K = Jypv'^dV, W = - Jy pr ■ V'^dV, 
M = (I/Stt) Jy B'^dV, S = -Jg[Pr + r- (prr)] ■ dS, and F = J^r -Tm ■ dS. The quantities 
J, U, K, W, S, M, F, and Tm denote, respectively, the moment of inertia, internal 
thermal energy, internal kinetic energy, gravitational energy, the sum of the thermal surface 
pressure and dynamical surface pressure, internal magnetic energy, the magnetic surface 
pressure, and the Maxwell stress-energy tensor. The internal thermal, kinetic, and magnetic 
energies are always positive. Other terms can be either positive or negative. In particular, 
the gravitational term can be positive in crowded enviro nments where the background 



Ballesteros-Paredes et al. 



20091 ). 



gravitational field dominates the core dynamics (see e.g.. 
In fact, as shown below, several cores have positive gravitational terms in our simulations, 
although no cores have positive gravitational terms for the initially magnetically subcritical 
case where the turbulent motions are not so strong because of the effects of magnetic cushion 
(see Fig. 16 of iNakamura fc Lill2008l ). The second term on the left-hand side of equation 



2]) denotes the time derivat ive of the flux of moment of inertia through the core boundar y. 



As is the standard practice (JTilley fc Pudritz 



2007 



Dib et al. 



2007 



Nakamura fc Li 



20081), 



we ignore the left-hand side of equation ([2]) in our discussion and consider a core to be in 
virial equilibrium if the sum U + K + W + S + M + F = 0. 

The equilibrium line is shown in Figure where the sum of the surface terms (5* + F) 
is plotted against the gravitational term (W), both normalized to the sum of the internal 
terms {U + K + M). For the surface term, the external kinetic term is generally much larger 
than the magnetic term (F). 

For all three models, the majority of the cores lie below the virial equilibrium line, 
and thus they are expected to be unstable to contraction. In addition, the surface term 
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appears more important than the gravitational term for most of the cores. In the absence 
of the magnetic field, about 10 — 20 % of the cores have the normalized surface term 
of \S + F\/\U + K + M\ > 5, whereas no cores have the normalized surface term of 
\S + F\/\U + K + M\ > 5 in the presence of the moderately-strong magnetic field. The 
magnetic field tends to reduce the contribution of the surface terms to the core dynamical 
state, although the surface term is still dominant even in the presence of the strong magnetic 
field. 

The core virial state is in contrast to the initially-magnetically subcritical case. For the 
initially-magnetically subcritical case, for most of the cores, the surface term [S + F) and 
gravitational term jW ) tend to be smaller than the sum of the internal terms, U + K + M 



(INakamura fc Li 



20081 ). and thus most of the cores are distributed around the virial 
equilibrium line. In addition, more massive cores tend to be distributed below the virial 
equilibrium line, and the importance of the self-gravity tends to depend on the core mass. 
The gravitational term tends to be more important for more massive cores for the initially 
magnetically subcritical case, implying that the gravitational collapse induced by ambipolar 
diffusion regulates the core evolution. In contrast, for the magnetically-supercritical cases 
studied in this paper, most of the cores are distributed further from the virial equilibrium 
line. The surface term (5* + F) tends to be larger than the sum of the internal terms, 
U + K + M for most of the cores. The core virial state shows no clear dependence on the 
mass. The surface term, which is typically dominated by the ambient supersonic turbulence, 
is still important even for massive cores. 



The above result is different from that of 



Tilley fc Pudritzl (120071 ). They found, in their 



ideal MHD simulations, a larger contribution from the surface term to the virial equation 
for a stronger magnetic field, which is oppose to the trend we found. This difference likely 
comes from the fact that they only followed the very early phase of cloud evolution, when 
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the initial supersonic turbulence has decayed significantly but star formation has yet to set 
in. Most of our cores are formed at much later times, after the decayed turbulence has been 
replenished by outflow feedback. Given this difference, it is perhaps not too surprising that 
the surface terms are much larger in our models than in theirs. This is particularly true for 
the weaker magnetic field case, because of a more rapid and violent star formation. 

In our simulations, most of the identified cores have the densities smaller than the 
critical density beyond which the Jeans condition is violated. For example, the Jeans 
condition is satisfied for 89 % (model Nl), 85 % (model Wl), and 95 % (model SI) of 
the cores presented in this subsection. For the stronger initial magnetic field, the mean 
density of the cores tend to be smaller due to the magnetic support, and therefore the 
Jeans condition is satisfied for a larger number of cores. We also performed a run with 
512^ for model SI, with the threshold density fixed, so that the Jeans condition is always 
satisfied in the entire computation box. To compare directly with the results presented in 
this subsection, we first resized the data obtained from the 512^ run to the 256^ grid data, 
and then applied the clumpfind method to the regridded data. We confirmed that the 
statistical properties of the identified cores are essentially the same as those presented in 
this subsection. For example, for all the cores identified from the 512^ data, the nonthermal 
3D velocity dispersions of the identified cores stay below 1 km s~^, and its mean is estimated 
to be 0.43 km s^^, comparable to that of the 256'^ data, 0.46 km s^^. 

3.4. Effects of Outflow Strength 

After exploring the effects of the magnetic field on the core properties, we now turn 
to the outfiow strength. In our model, the strength of the outfiow feedback is specified 
by the dimensionless parameter, /. This parameter is highly uncertain because it is 
generally difficult to determine the physical properties of molecular outflows accurately 



from observationa 
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data. In this paper, we adopt / = 0.4 as a fiducial value, following 



Matzner Sz McKed (|2000| ). Recently, Nakamura et al. (2011) discovered many molecular 
outfiow lobes toward an extremely-young cluster-forming clump. Serpens South, on the 
basis of CO ( J = 3 — 2) observations, and derived the physical properties of the identified 
molecular outfiows. They roughly estimated / ~ 0.1 — 0.3 for these outfiows, assuming that 
the outfiow gas is optically-thin. If a substantial amount of the outfiow gas is optically-thick, 
the outfiow parameter / would be larger. In the following, we compare the core properties 
of three models with different outflow strengths: / = 0, 0.1, and 0.4. The flducial value for 
plasma (3 of 0.2 is adopted for these models. 

In Figs. [6^ through int, we compare the velocity dispersion-radius relations for the 
moderately-strong magnetic models with no outflow feedback (/ = 0), weak outflow 
(/ = 0.1), and strong outflow (/ = 0.4). In the model with no outflow feedback (model 
S3), the number of identifled cores is smaller than in the other two models, and the mean 
core mass tends to be slightly larger. This difference may be due to the fact that the 
ambient turbulent motions are weaker without any outflow feedback and thus the turbulent 
fragmentation does not form smaller-scale structures as efficiently. Figure |6] indicates that 
the velocity dispersion-radius relation does not depend strongly on the strength of the 
outflow feedback. For all three models, there is no clear relationship between the velocity 
dispersion and the core radius. Thus, it appears difficult to discriminate the effects of 
protostellar outflows based on the velocity dispersion-radius relation alone. 

The surface term-gravitational term plots in Figs. [7^ through [Tt indicate that the 
dynamical states of the cores depend somewhat on the outflow strength. In the absence of 
the outflow feedback, more cores are located near the dashed line (where the two terms are 
equal). The gravitational term appears more important for this model than for the other 
two, although the surface term is clearly still important. This is due, at least in part, to 



-19- 

the fact that the global gravitational infall can generate fast motions even in the absence of 
protostellar outflow driven turbulence. 

It is difficult, however, to gauge the effects of the outflow feedback from the virial 
parameters of the cores, which are often used in observational studies of the molecular cloud 
cores in star-forming molecular clouds. The difficulty is illustrated in Figs. ^ through 
[Ht, where we show the virial parameter-mass relations for the three models with different 
outflow strengths. In all cases, the distributions of the virial parameters are similar, 
following a power-law of a^h oc Mc with rather small dispersions. 

3.5. Outflow Direction 

In this subsection, we investigate how well the outflow axes are aligned with the global 
or initial magnetic fleld direction in parsec-scale cluster forming clumps. In our simulations, 
the outflow direction is set to the local magnetic fleld direction at the position of each 
formed star. This should be a reasonable approximation since the outflows are most likely 



driven along the 



Matsumoto et al. 



■ocal magnetic fleld by fleld line twisting due to disk rotation, (e.g.. 



20061). 



Figure [9] shows the histogram of the outflow direction as a function of the angle 
between the outflow axis and the initial magnetic fleld direction, for the model with a 
moderately-strong magnetic fleld (model SI). All the outflows created by the time when 
SFE ^ 0.16 are shown. The outflow direction shows a broad distribution, taking its 
maximum at the mean angle of about 30°. The conclusion is that, although there is some 
preference for the outflows to align with the global magnetic fleld direction, the alignment 
is not as strong as one may naively expect. This is because most of the dense cores are 
formed through outflow-induced turbulent compression, which can push the magnetic fleld 
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on the core scale away from the global field direction. On smaller scales, the misalignment 
in field direction can be further amplified by gravitational collapse. The lack of a complete 
alignment between the outfiow axes and the global field direction is thus expected. It can 
not be used to argue against the presence of a moderately-strong global magnetic field. 

4. Discussion 



4.1. Comparison with Core Observations 

Here, we compare the core properties obtained in Section [3] with observations. We use 
the core data toward two nearby parsec-scale cluster-forming clumps: the p Ophiu chus Main 



Cloud and NGC 1333 in Perseus. In the p Ophiuchus region. 



Maruta et al. 



mm identified 



68 H^^CO+ (J = 1 — 0) cores by applying clumpfind to the 3D posit ion-position- velocit y 



Walsh et al 



mm 



cube data obtained using Nobeyama 45 m telescope. In NGC 1333, 
identified 93 N2II+ ( J = 1 — 0) dense cores, again using clumpfind. Since both the H^^CO"'" 
( J = 1 — 0) and N2H"*" (J = 1 — 0) transitions have critical densities of order ~ 10^ cm^'^, the 
H^^CO"'" and N2H+ cores should have densities comparable to our simulated cores, which 
were identified using a threshold density ~ 10^ cm~'^. Since the moderately-strong field 
model yields a core velocity dispersion that is in better general agreement with observations 
than the non-magnetic or weak-field model, we will concentrate on comparing its cores with 
the p Oph and NGC 1333 data. 

In Fig. [TUh . we c ompare the ID FWH M nonthermal velocity width-radius relations 



for the p Oph cores of 



Maruta et al. 



(J2OIOI ) and our simulated cores (see Section [3]). The 
3D velocity dispersions discussed in Section [3] was converted into the ID FWHM velocity 
widths using c^Vnt = (81n2/3)^/^AVNT- Our simulated cores tend to be somewhat smaller 
than the p Oph cores; the average radius of the p Oph cores is estimated at 0.04 pc, about 
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1.7 times that of model SI. Aside from the size difference, the observed and simulated 
data are almost indistinguishable. Neither shows a clear dependence of linewidth with 
size, except for a weak trend for the linewidth to increase slightly with radius. Neither 
follows the Larson's linewidth-size relation. The observations are also consistent with the 
cores formed from initially-magnetically subcritical clouds (see Fig. 18 of 



Nakamura fc Li 



1 of 



2008). We did not plot the NGC 1333 cores on Fig. [TUh because the core radius in Table 



Walsh et al 



m 



( 120071 ) was quoted in single significant digit, which artificially introduces 
a large scatter in the plot. Nevertheless, inspection by eye reveals no clear correlation 
between the velocity width and core radius, broadly consistent with the data for both the p 
Oph cores and the simulated cores of model SI. 

Figure [TOb compares the virial parameter-mass relations of the simulated and obse rved 



cores. The red an d blue dots, and b lack crosses denotes the p Oph (JMaruta et al. 



NGC 1333 cores flWalsh et al 



2010h and 



20071 ). and the cores of model SI, respectively. For the NGC 



1333 cores, we recalculated the virial masses by including the thermal contribution 
the gas with a temperatu re of 15 K and adop ting a distance of 235 pc f.Hirota et al. 



instead of 300 pc used in 



Walsh et al. 



rom 



20081), 



(120071 ). We also recalculated the LTE masses by 



adopting the N2H"'" fractional abund ance relative to Ho of Xn2h+ ~ 3 x 10 ^° derived from 



observations toward p Oph because 



Walsh et al. 



(I2OO7I ) determined their adopted N2H+ 



fractional abundance such that the points are distributed evenly about the line of equality 
(Mlte = Mvir) in their MuTE-M^h plot. This N2H+ fract ional abundance is deter mined 



by taking aver age of the values ob tained toward Oph A ( JDi Francesco et al. 



B2 sub clumps (JFriesen et al. 



2OO4J ) and 



2010l ). Figure [TUb indicates that the virial parameter-core 



mass relation of the nearby cluster-forming clumps follows roughly a single power-law of 

—2/3 

Ovir oc Mc with a small dispersion, which is again in good agreement with model SI that 
has a moderately-strong magnetic field. A dynamically important magnetic field may have 
played a role in core formation in these regions. 
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Besides the p Oph region, 



Maruta et al. 



(120 lOl ) found that the velocity dispersions of 



H^^CO"'' cores are almost independent of the core radius in Orion A as well, even though the 
latter is more than a factor of 3 farther away. The coarser spatial resolution for Orion means 
that its cores may contain (smaller, p Oph-like) cores blending together. If this is true, 
then the nearly flat velocity-width-radius relation observed in both regions suggests that 
the inter-core motions among the neighboring cores are almost cor nparable to the inte rnal 



Andre et al. 



( 120071 ) who 



motions in the individual cores. Such a feature was pointed out by 
measured the velocity difference among neighboring cores using N2H+ (J = 1 — 0). If the 
velocity-width-radius relation is flat (AV oc constant) and the core mass is proportional 
to -Rj?, then the virial param eter is scaled as a^w oc R rAV'^ IJGMr) oc Mc , a similar 



power-law to that derived by 



Bertoldi fc McKee 



( I1992I ). This again suggests the importance 



of ambient turbulent pressure in dynamics of the cores as discussed in Section El 



4.2. Amplification of Magnetic Field 

To isolate the effect of magnetic field on clustered star formation, we present in Fig. 
ITTk the time evolution of star formation efficiency (SFE) for three models with different field 
strengths but no outfiow feedback. Our simulations show that, as expected, the magnetic 
field tends to slow down star formation, especially at relatively late times. The rate of star 
formation remains rather high, however, particularly at early times. For example, from 
the formation of the first star to the time when the SFE reaches ~ 16%, the average star 
formation efficiency per global free-fall time (SFRg; Krumholz & McKee 2005) is estimated 
to be 29 %, 27 %, and 21 %, respectively, for models N2, W2, and S2. These valu es are 



larger than the observationally inferred values of a few % ( jKrumholz fc Tanl 120071 ). The 
large SFRg in the absence of outfiow feedback is in agreement with previous studies (see 
also the Appendix for the analytic formula of SFRg). For example, using MHD SPH 
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simulations, iPrice fc Batd ( 120081 ) and iPrice fc Batd ( 120091 ) followed the evolution of 50 M© 
clumps with 0.25 < /3 < oo until t < 1.5tff^ci- Although their initial model clump mass is too 
small compared to nearby cluster-forming clumps like p Oph and NGC 1333 (which have 
masses of order 10^ M©), their star formation efficiencies per free-fall time are estimated 
to be over 10 %. These simulations imply that it is difficult to reduce the star formation 
rate to the observed level by a moderately-strong magnetic field alone, and other factors 
are needed to significantly retard star formation in cluster-forming clumps. 



One way to slow down star formation is the feedback from forming stars. 



Price fc Bate 



(120091 ) found that radiative feedback from forming stars does not change the global star 
formation efficiency in a parent clump much, although it significantly suppresses the 
small-scale fragmentation by increasing the temperature in the high-density material near 
the protostars. Figure [TTb shows that the star formation rate is greatly reduced by the 
inclusion of outflow feedback, eve n in the presence of a relatively weak magnetic field, in 



agree ment with previous studies (JLi fc Nakamura 



2006 



Nakamura fc Li 



20071 : 



Wang et al. 



2010l ). The average star formation efficiency per global free-fall time is estimated to be 8 %, 
6 %, and 4 % for model Nl (no magnetic field), Wl (weak field), and SI (moderately-strong 
field), respectively. For the model with the moderately-strong field, the SFRfj is more or 
less comparable to the observed value. 

The reason for the large reduction in SFRfj can be seen in Fig. [121 where the magnetic 
energy is plotted against the evolution time, for models Wl and SI. As shown in Fig. 
[T^ . for the moderately-strong field case, the total magnetic energy is dominated by 
the background uniform field that does not contribute to the force balance in the initial 
cloud. We therefore illustrate in Fig. [T2b the time evolution of the magnetic energy 
stored in the distorted component that was amplified by supersonic turbulence. Here, we 
computed the magnetic energy stored in the distorted component by subtracting the initial 
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magnetic energy from the total magnetic energy. For comparison, we also plotted in Fig. 
[T2b the evolution of the total kinetic energy of the clump. Since the kinetic energy is 
often dominated by the unbound high-velocity gas associated with the active outflows, we 
approximate the total kinetic energy of the clump as the kinetic energy of the gas whose 
velocity is smaller than 10 Cg. For both models, the amplified components tend to increase 
with time. For the model with the moderately-strong field, the amplified component begins 
to oscillate around a level value after a free-fall time. For both models, the magnetic energy 
of the amplified component becomes comparable to the kinetic energy of the dense gas by 
the end of the compu tation, indicating that a quasi-equipartition has been reached (see also 



Federrath et al. 



2OIII ). For the weak field case, the amplified component is more important 
than the initial uniform field, resulting in a significantly-distorted magnetic field structure. 
Such a highly-distorted magnetic field can be seen in the 3D bird's eye view of the density 
and magnetic field distributions of model Wl, which is presented in the left panel of Fig. 
In contrast, the global magnetic field is well-ordered for the case of moderately-strong 



initial magnetic field, which is shown in the right panel of Fig. [131 



4.3. Polarization Maps 



Polarization maps of submillime ter thermal dust emission have recent l y been obtained 
for n earby star forming regions (e.g 



Houde et al 



2004J : 



Girart et al 



20091 



Matthews et al. 



20091 ) . Here, we present the polarization maps derived from the simulation data for the two 
magnetically-supercritical models with different initial magnetic field strengths (models Wl 
and SI which have the magnetic flux-to-mass ratios of A = 4.3 and 1.4, respe ctively). We 



Padoan et al. 



compu ted the polarized thermal dust emission from the MHD model following 

(I2OOII ). We neglect the effect of self-absorption and scattering because we are interested in 

the thermal dust emission at submillimeter wavelengths. We further assume that the grain 
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properties are constant and the temperature is uniform. The polarization degree is set such 
that the maximum is equal to 15 %. Figure UM shows the dust polarization maps calculated 
from models Wl and SI. 

Only a small portion of the computation box is shown in each panel of Fig. [TH As 
expected from Fig. [T31 in the presence of a weak magnetic field, the spatial distribution 
of the polarization vectors has relatively large fluctuations. The polarization degree tends 
to be smaller in several parts where the magnetic fields are strongly distorted. In the 
densest parts, the polarization vectors appear more or less parallel to the local elongated 
structures or dense filaments, whereas there is no clear correlation between the vectors 
and density distribution in less dense parts. This may be due to the fact that the local 
dense filaments are created by turbulent compression that preferentially enhances the 
magnetic field component transverse to the shock plane. In contrast, in the presence of 
the strong magnetic field, the filamentary structure is prominent and the filament axes 
tend to be perpendicular to the polarization vectors that are almost parallel to the initial 
magnetic field direction, suggesting that the gas flow is channeled preferentially by the 
strong magnetic field to form the filaments. The polarization observations can thus provide 
a handle on the magnetic field strength of cluster-forming clumps. 



Recent polarization observations of cluster-forming clumps s how that the 



magnetic fie 



Houde et al. 



d lin e s are more or less spatia 



2004 



Sugitani et al 



2010 



ly well-ordered (e.g., 



201ll ). For example. 



Davis et al 



Sugitani et al 



glob al 



2000 



teoioh found 



that the Serpens cloud core is threaded by a hou r-glass shaped well-or dered magnetic 



Sugitani et al. 



(1201 ih found that the 



field and is elongated in the cross-field direction. 

Serpens South filamentary infrared dark cloud appears to be threaded by more or less 

straight global magnetic field. The observed spatially-ordered magnetic field structures 



imply that the magnetic fields in the nearby cluster-forming clumps are likely to be at least 



-26- 
moderately-strong and djTiamically important. 

4.4. Core Mass Function 

In Figs. [T5k . [T5b . and [T5b . we plot the core mass functions (CMFs) for our identified 
cores for models Nl, Wl, and SI, respectively. To enlarge the number of sample cores, 
we added up all the cores identified at the three stages when SFE has reached 0.08, 0.12, 
and 0.16. In addition, we added up the cores identified from the runs using the different 
initial turbulent realization for which the different random numbers were used for both the 
amplitude and phase. The total number of cores so identified is over 500 for each model. 

The overall shapes of the CMFs are similar in all three cases, suggesting the shape 
of the CMF is relatively insensitive to the initial magn etic field strength. T his is different 
from the turbulent fragmentation scenario proposed by 



Padoan et al. 



( 1200 if ) where a weak 



magnetic field is required to reproduce the CMF that is similar in shape to the Salpeter 
IMF (see also lLi et al.il2010i ). In our simulations, a moderately-strong magnetic field can 
produce a CMF that resembles the Salpeter IMF as well. 

We note two interesting features of the computed CMFs. First, there appear to be 
a lack of massive cores compared to the Salpeter stellar IMF. This is not necessarily a 
problem, because it is unclear whether the nearby cluster-forming clumps such as p Oph 
and NGC 1333 that we aim to simulate would ever produce massive stars. Alternatively, 
massive stars can grow from initially less massive cores, fed by global collapsing flows 
towards the bottom of the clump potential (Smith et al. 2008; Wang et al. 2010). Second, 
the turnover at the low mass end is less clear in the CMFs than in the Kroupa or Chabrier 
IMF. This may be because, at the low mass end, most of the cores are gravitationally 
unbound, and may not form stars. Further studies are required to determine the relation 
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between the CMF and IMF. 

5. Summary 

We have performed a set of 3D MHD simulations of cluster formation taking into 
account the effects of protostellar outflow feedback, and identified dense cores by applying 
a clumpfind algorithm to the simulated 3D density data cubes. The main results are as 
follows. 

1. Dense cores do not follow Larson's linewidth-size relation. We find that the velocity 
dispersions of dense cores show little correlation with core radius, irrespective of the 
strength of the magnetic field and outflow feedback. In the absence of a magnetic field, 
the majority of the cores have supersonic velocity dispersions, whereas in the presence of a 
moderately-strong magnetic field, the cores tend to be subsonic or at most transonic. 

2. We find that most of the cores are out of virial equilibrium, with the external pressure 
due to ambient turbulence dominating the self-gravity. The core formation and evolution is 
largely controlled by the dynamical compression due to outflow-driven turbulence. Such a 
situation is contrast to the strongly-magnetized (magnetically subcritical) case, where the 
self-gravity plays a more important role in the core dynamics, particularly for massive cores. 

3. Even an initially- weak magnetic field can retard star formation significantly, 
because the field is amplified by supersonic turbulence to an equipartition strength. In 
such an initially weak field, the distorted field component dominates the uniform one. In 
contrast, for a moderately-strong field, the uniform component remains dominant. Such 
a difference in the magnetic structure can be observed in simulated polarization maps 
of dust thermal emission. Recent polarization measurements show that the field lines in 
nearby cluster-forming clumps are spatially well-ordered, indicative of a moderately-strong. 
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dynamically-iinportant, field. 
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A. Star Formation Rate in A Cluster Forming Clump 



Using 3D MHD n umerical simulations of cluster formation, iLi &: Nakamural (120061 ) and 



Nakamura &: Lil (120071 ) demonstrated that the protostellar outflow-driven turbulence can 
keep a pc-scale, cluster-forming clump close to a virial equilibrium long after the initial 
turbulence has decayed away. Here, we derive an analytic formula of star formation rate in a 
cluster-forming clump that keeps its virial equihbrium by the protostellar outflow feedback. 

Numerical simulations of protostellar turbulence indicate that the dissipation rate 
of the turbulence momentum, dPturh/dt, balances the momentum injection rate by the 
protostellar outflow feedback, dPout/dt, so that the clump can be kept close to a virial 
equilibrium, 

(Al) 



turb 



out 



dt 



dt 



The dissipation rate of the turbulence momentum, dPturh/dt, can be written as 



dR 



turb 



M,,K 



dt 



a- 



cl >^vir 



where a ~ 0.5 (iMac Low 



^dif 



19991 ). The virial velocity, Kir, is given by 

V\r\r — 



3 „ GMei 



(A2) 



(A3) 



The momentum dissipation time is tdiss = Rc\/Vvh- The dimensionless parameter of order 
unity, a, measures the effects of a nonuniform or nonspherical mass distribution. For a 



uniform sphere and a centrally-condense sphere with p oc r 



1 and 5/3, respectively. 



In the following analysis, we adopt a = 5/3 because the cluster-forming clumps tend to be 
centrally-condensed. We take into account the magnetic support by multiplying the virial 
velocity by a factor /b, where < /b < 1. 



The momentum injection rate by the protostellar outflows is given by 

"-Tout 



dt 



^SFR X /wKv 



(A4) 
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where V^ is the flow speed and esFR is the star formation rate. Using equations flAl|) and 
flA4p . the star formation rate is rewritten as 



esFR = -aafsU K —52- 



3 

5^ 



(A5) 



If we normahze espR to Md/tg, then the star formation efficiency per free-fall time is given 

by 



SFR, = esPK^ - 0.0125a (A^ (^) (lO^I^) (^) V 5 x 10-cm- 



^cl 



1/2 



(A6) 



where the surface density is S = Md/nRli and the free-fall time is tg = ^y3n/32Gp. 



Figure [16] shows the dependence of SFRg on the mass and radius, obtained from the 
above equation. The crosses an d diamonds i ndicat e the cluster-forming clum ps observed 



in the C^^O ( J = 1 - 0) fine by 



Ridge et al. 



(12003h and 



Higuchi et al. 



( l2009l ). respectively. 



Our model suggests that the star formation rate per free-fall time ranges from 1 % to 5 % 
for the observed cluster-forming clumps in the solar neighborhood when the protostellar 
outflow feedback maintains the supersonic turbulence in the clumps, indicating that it takes 
about (2 — 10) iff for the star formation efficiency to reach about (10 — 20) %. 
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Table 1. Model Characteristics 



Model /3 /=" Remark 

Nl oo 0.4 no B field 

N2 oo 0.0 no B field, no feedback 

Wl 2 0.4 weak B field 

W2 2 0.0 weak B field, no feedback 

81 0.2 0.4 moderately-strong B field 

52 0.2 0.0 moderately-strong B field, no feedback 

53 0.2 0.1 moderately-strong B field, weak feedback 



^outflow strength 
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Table 2. Summary of the Physical Properties of the identified cores 

Property Minimum Maximum Mean'^ Median 

model Nl 

Re (pc) 0.013 0.041 0.020 ± 0.005 0.020 

Mc (Mq) 0.075 3.51 0.67 ± 0.67 0.40 

Al^NT (km s-^) 0.13 4.53 1.06 ± 0.77 0.86 

model Wl 

i?c (pc) 0.013 0.034 0.020 ± 0.006 0.019 

Mc (Mq) 0.076 3.90 0.62 ± 0.77 0.31 

AKnt (km s-i) 0.13 14.4 0.99 ± 1.66 0.47 

model SI 

i?e (pc) 0.013 0.048 0.024 ± 0.008 0.023 

Mc (Mq) 0.073 4.39 0.86 ± 1.03 0.47 

AV^NT (km s-i) 0.086 2.64 0.46 ± 0.52 0.31 

^With standard deviation 

Note. — The gas temperature is assumed to be T = 20 K, corresponding to the sound 
speed of 0.266 km s^^ 
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Fig. 1. — Column density snapshots on the x-z plane for (a) model Nl (t = 4.6tff^c ~ 0.9 
Myr), (b) Wl {t = 6.5tff,c ~ 1.2 Myr), and (c) SI (t = 7Ats,c ~ 1-4 Myr) at the stage when 
16 % of the total mass has been converted into stars. The positions of the formed stars are 
overlaid with the dots. The units of length are the central Jeans length Lj = (vrc^/Gpo)^^^ — 
0.17(T/20K)^/^(r7,H2,o/2.69 x lO'^cm^^)^^/^ pc. The initial magnetic field direction is parallel 
to the X-axis. In panels (d) through (f), the positions of the identified cores are overlaid with 
the triangles on the same column density images as those shown in panels (a) through (c). 
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Fig. 2. — Histograms of some physical quantities of the identified cores. They are plotted at 
the same stages as those of Fig. [TJ Panels (a), (b), and (c) are the histograms of the core 
radius for models Nl, Wl, and SI, respectively. Panels (d), (e), and (f) are the histograms 
of the core mass for models Nl, Wl, and SI, respectively. Panels (g), (h), and (i) are the 
histograms of the 3D velocity dispersion of the nonthermal component for models Nl, Wl, 
and SI, respectively. The vertical dashed lines indicate the sound speed with gas temperature 
of 20 K. 
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Fig. 3. — Nonthermal 3D velocity dispersions as a function of core radii for three models 
with different magnetic field strengths: (a) model Nl, (b) model Wl, and (c) model SI. 
They are plotted at the same stages as those of Fig. [TJ The sound speed is Cg = 0.266 km 
s"\ corresponding to T = 20 K. 
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Fig. 4. — Virial parameters as a function of core masses for three models with different 
magnetic field strengths: (a) model Nl, (b) model Wl, and (c) model SI. They are plotted 
at the same stages as those of Fig. [31 
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Fig. 5. — Relationship between the sum of the two surface terms, S+F, and the gravitational 
term, W , in the virial equation. They are normalized to the sum of the internal terms, 
U + K + M. The sohd line indicates the virial equilibrium, U + K + W + S + M + F = Q. 
They are plotted at the same stages as those of Fig. [31 The dashed line indicates the line 
below which the surface term is larger than the gravitational term. For the cores that lie 
below the solid line, the left-hand side of the virial eq. [2] is negative and thus expected to 
be bound. All others that lie above the solid line are unbound and expected to disperse 
away, if they do not gain more mass through accretion and/or merging with other cores, 
or reduce internal support through turbulence dissipation. A red cross, black triangle, and 
blue square indicates a core with mass of Mc < 0.75 (Mc), 0.75 (Mc) < Mc < 1.25 (Mc), and 
Mc > 1.25 (Mc), respectively. The close-up of Fig. [St is presented in Fig. ^. 
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Fig. 6. — Same as Fig. |3]but for models (a) S3, (b) S2, and (c) SI. 




-2-10 1 2 

W/(U+K+M) 




-2-10 1 2 

W/(U+K+M) 




-2-10 1 

W/(U+K+M) 



Fig. 7. — Same as Fig. |5]but for models (a) S3, (b) S2, and (c) SI. 
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Fig. 8. — Same as Fig. Hlbut for models (a) S3, (b) S2, and (c) SI. 
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Fig. 9. — Histogram of the outflow direction relative to the initial magnetic field direction 
for the model with a moderately-strong magnetic field (model SI). 
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Fig. 10. — (a) the relationship between the ID FWHM nonthermal velocity width and core 
radius for the dense cores in the p Ophiuchi Main Cloud and our cores identified from the 
moderately-strong magnetic field model (model SI). The red dots and black crosses denote 
the p Op h cores and ou r cores of model SI, respectively. For the p Oph cores, we used the 



results of 



Maruta et al 



(120 inh who identified the cores using H^^CO+ (J = 1-0) data. 
For our cores of model SI, we recalculated the ID FWHM velocity widths from the 3D core 
velocity dispersions derived in Section El (b) The virial parameter-core mass relation. The 
red dots and black crosses are the same as those of panel (a). The blue open circles d enote 



the cores of NGC 1333 identified using N2H+ ( J = 1 - 0) emission flWalsh et al.lboOTI ). For 



p Oph cores and NGC 1333 cores, the LTE masses are adopted as the core masses. For 



20071) by 



the NGC 1333 cores, we recalculated the LTE masses derived by ( iWalsh et al. 
assuming the average N2H"'' fractional abundance relative to H2 obtained toward p Oph of 
-^N2H+ ~ 3 X 10^^° and distance of 250 pc. See the text for detail. 
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Fig. 11. — (a) Time evolution of star formation efficiencies for three models with no outflow 
feedback. The three models have different initial magnetic field strengths: models N2 (/3 = 
oo), W2 (/3 = 2), and S2 (/3 = 0.2). The evolution time is normalized to the global clump 
free-fall time, tfj^ci = 0.49(L/1.5pc)(20K/T)^'^ Myr. (b) Same as panel (a) but for models 
with outflow feedback (models Nl, Wl, and SI). 
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Fig. 12. — (a) Time evolution of total magnetic energy for two models with different initial 
magnetic field strengths. Thick and thin lines are for the models with moderately-strong 
(model SI) and weak magnetic field (model Wl), respectively. The evolution time is nor- 
malized to the global clump free-fall time, tg^ci = 0.49(L/1.5pc)(20K/T)^/^ Myr. (b) Time 
evolution of magnetic energy of amplified component (solid lines) and kinetic energy of tur- 
bulence (dotted lines). The magnetic energy of the amplified component tends to approach 
the kinetic energy for each model. 
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Fig. 13. — 3D view of the density and magnetic field distributions for (a) tlie weakly- 
magnetized model (model Wl) and (b) the strongly-magnetized model (model SI) at the 
stage where the star formation efficiency has reached 0.16. The color images indicate the 
isodensity surfaces. The white curves are the magnetic field lines. For both the panels, only 
the central L/2 x L/2 x L/2 region is shown. 
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Fig. 14. — (a) Polarization maps of (a) model Wl (weak magnetic field) and (b) model 
SI (moderately-strong magnetic field). The length of a polarization vector is proportional 
to the degree of polarization, with the longest vector corresponding to P = 15%. Only one 
polarization vector is plotted for every four computational cells. The color contour shows the 
column density distribution. The initial magnetic field lines are parallel and perpendicular 
to the horizontal line for panels (a) and (b), respectively. The units of length are the central 



Jeans length Lj = {ncl/GpoY/^ ~ 0.17(T/20K)i/2(r2jj^ o/2.69 x lO^cm" 



-3N-1/2 



pc. 
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Fig. 15. — Core mass functions for (a) model Nl, (b) model Wl and (c) model SI. A 
power-law of AN/ AM oc M~^ is plotted for comparison. 
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Fig. 16. — Contours of star formation rate per free-fall time predicted by the outflow- 
regulated cluster formation model on the mass-radius diagram. The dashed contours are 
labeled by values of SFRg. The dotted contours indicate the constant column density at 
10^^ cm~^ [bottom line), 10^^ [intermediate line), and 10^^ cm~^ ( upyer line). The cr osses 



and triangles i ndicat e the cluster-forming dense clumps identified by 



Higuchi et al 



Ridge et al. 



fl2003[ ) and 



( 120091 ) in CO (J = 1 — 0), resp ectively. The ast e risks indicate the dense 



clumps in the infrared dark clouds identified by 



Rathborne et al. 



(120031). For most of the 



dense clumps, the predicted SFRg ranges from 1 % to 5 %, indicative of slow star formation. 



